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Abstract 

The elementary vortex pinning potential is studied in unconventional superconduc- 
tors within the framework of the quasiclassical theory of superconductivity. Numeri- 
cal results are presented for d-, anisotropic s-, and isotropic s-wave superconductors 
to show explicitly that in unconventional superconductors the vortex pinning poten- 
tial is determined mainly by the loss of the condensation energy in bulk due to the 
presence of the pinning center, i.e., by the breakdown of the Anderson's theorem. 
It is found that the vortex pinning energy in the d-wave pairing case is 4 - 13 times 
larger than those in the s-wave pairing cases. This means that an enhancement of 
pinning effect in unconventional superconductors occurs due to the breakdown of 
the Anderson's theorem. The case of a chiral p-wave superconductor is also inves- 
tigated in terms of the vortex core states subject to the Andreev reflection, where 
important is whether the vorticity and chirality are parallel or antiparallel. 
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1 Introduction 



Much attention has been focused on the vortex pinning in type-II superconduc- 
tors. The vortex pinning in the superconductors under magnetic fields plays 
an important role on various vortex-related quantities and phenomena such 
as the critical current, hysteresis of the magnetization, and quantum vortex 
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tunneling. The problem of the vortex pinning is categorized in two aspects. 
One is the elementary vortex pinning force, which is the interaction between a 
vortex and a single defect. The other is the summation problem, which is how 
the elementary vortex pinning forces add up to a holding force on the elastic 
vortex lattice under a distribution of defects [1]. A study of the elementary 
vortex pinning is the step necessary for the understanding of the phenomena 
related to the vortex pinning in the superconductors under magnetic fields. 

Traditionally, the mechanism of the elementary vortex pinning is considered 
to be that a defect locally prohibits superconducting condensation so that it 
attracts a normal region of a vortex core in order to avoid the loss of the 
condensation energy. Thuneberg et al. [2,3] advanced the understanding of 
the mechanism of the elementary vortex pinning, taking account of a nonlocal 
effect that a defect as a scattering center scatters the quasiparticles. This 
mechanism of the vortex pinning is that the quasiparticle scattering by the 
defect inside a vortex core helps to avoid suppressions of the order parameter 
around the vortex center up to distances of the order of the coherence length. 

In this paper, we investigate the elementary vortex pinning potential in un- 
conventional superconductors, following the analysis started by Thuneberg et 
al. [2,3] for conventional s-wave superconductors. The unconventional super- 
conductivity is recently proposed for many superconductors such as high-T^ 
cuprates, organic conductors, and heavy-fermion compounds. For example, 
the d-wave superconductivity is believed to be realized in most of the high- 
Tc cuprates, and the chiral p-wave one is expected in a non-copper-layered 
perovskite superconductor Sr2Ru04 [4] (although it seems to be open to fur- 
ther discussion [5]). We calculate numerically the elementary vortex pinning 
potential for d-, anisotropic s-, and isotropic s-wave pairings, and compare 
them each other to clarify what could be one of the characteristic properties 
specific to the superconductors with the unconventional order parameter. We 
also present an analytical consideration of a novel chirality-dependent vortex 
pinning in a chiral p-wave superconductor. 



2 Formulation 

To investigate the vortex pinning, we use the quasiclassical theory of super- 
conductivity [6] . We start with the Eilenberger equation for the quasiclassical 
Green function in the absence of the pinning. 
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iVpk ■ V^imt + [iUJnh - A, Qi^t] = 0, (2) 

where fa is the Pauh matrix and the order parameter 



A(r.k) = | Jlr. ^<'„-^'). (3) 



-A*(r,k) 

The Eilenberger equation (2) is supplemented by the normahzation condition 
^imt(^<^n) r, k)^ = —TT^l, and the commutator is [a,b] — ah — ha. The symbol 
r = (r cos 0, r sin 0) denotes the center of mass coordinate of the Cooper pairs 
and k = (cos ^, sin ^) denotes the relative coordinate of them. The cylindrical 
Fermi surface is assumed. We use units in which % = = ^■ 

Following Thuneberg et al. [2,3,7], the effect of the pinning is introduced as 
follows to the quasiclassical theory of superconductivity. The quasiclassical 
Green function g in the presence of a point-like non-magnetic defect situated 
at r = R is obtained from the Eilenberger equation 



ivY^s. ■ Vg + [iunh - A,g]^ [i, ^imt](^(r'), (4) 
and the t matrix due to the defect 



t{iu;n, r') = — [1 Nov{gimt(iuJn, r', k))ej , (5) 

where r' = r-R, the denominator D = l + (7riVot;)^[(5'imt)e + (/imt)e(/iLt)e], the 
symbol (■ ■ ■)6» = / • ■ ■d6/27i, the normal-state density of states on the Fermi 
surface is Nq, and we assume the s-wave scattering v to obtain Eq. (5). The 
cross section of the defect is given by atr — {Air/kp) sin^ 5o with the scattering 
phase shift Sq, and tan^o = —ttNqv. 

The free energy in the presence of the defect is, at the temperature T, given 
as [2,3,7] 



1 

5n{K) = NoT j dXY, j dk j dvTi[5gx^bl (6) 



where 5g\ = g — gimt is evaluated at A = AA5, and A5 is the order parameter 
in the absence of the defect. Equation (6) represents the difference in the free 
energy between the states with and without the defect, and then gives the 
vortex pinning potential 5i7(R). 
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Fig. 1. The vortex pinning potential as a function of the distance R between the 
vortex center and the defect. The defect is situated on the x axis (in the direction 
of = 0). 



3 Results 



3.1 d-, anisotropic and isotropic s -wave pairings 



We calculate the vortex pinning potential (5r2(R) numerically with the or- 
der parameter around a vortex, A(,(r, k) = F(6')Abcs(^) tanh[r/^(T)] exp(i0), 
(note that the vortex center is situated at r = 0). Here, the k-dependent func- 
tion F{9) — 1 for isotropic s-wave, -\/2|cos2^| for anisotropic s-wave, and 
-s/2cos2^ for d-wave pairings. Abcs(^) is the BCS order parameter amplitude 
and e(T) the coherence length. {^{T) = • [Ao/ABcs(r)], Aq = Abcs(O), and 

^0 = Vf/Aq.) 

In Fig. 1, we show the numerical results for the Born hmit {Sq 0) and the 
unitary limit {So = vr/2). The defect is situated on the x axis (in the direction of 
= 0). In the case of the isotropic s-wave pairing, when the defect is situated 
in bulk, i.e., far away from the vortex center {R = |R| oo), the difference in 
the free energy between the states with and without the defect, Sfl{R), is equal 
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to zero. This is because the Anderson's theorem [8] is vahd for homogeneous 
systems in isotropic s-wave superconductors and the non-magnetic defect has 
no influence on the free energy. On the other hand, in the anisotropic pairing 
cases, 6n{R — > oo) is finite and positive as seen in Fig. 1, which is the loss 
of the condensation energy in bulk due to the presence of the defect and to 
the breakdown of the Anderson's theorem. The condensation energy loss in 
bulk 6Q{R oo) contributes to the depth of the vortex pinning potential 
6Q{R), i.e., the vortex pinning energy. While the contribution of 6Q{R — > oo) 
to the vortex pinning energy was first pointed out by Thuneberg et al. [7] 
in the context of an ion in p-wave superfiuid ^He-B, we performed here a 
numerical calculation explicitly for the vortex pinning potential in the rf-wave 
superconductors. It is noted from Figs. 1(c) and 1(d) that, at T = 0.8Tc (high 
temperature), the vortex pinning energy in the d-wave case is dominantly 
determined by the condensation energy loss in bulk 5Q.{R — > oo), and is 4 - 
13 times larger than those in the s-wave cases. 

3.2 chiral p-wave pairing 

For the chiral p-wave pairing d = z(kx ± iky) [4], we analytically investigate 
the vortex pinning potential. In the case of that chiral p-wave pairing, there 
exist the following two possible forms for the order parameter A{r,(f);6) (= 
A(r, k)) of axisymmetric vortex [9-11]: (i) A^^(r, 0;6') = A_|_(r) exp[i(^ — 
0)]-|-A_(r) exp[i(— ^-1-0)], where the chirality and vorticity are antiparallel; (ii) 
A^+(r,0;6l) = A+(r) exp[i(6' + c6)] + A_(r) exp[i(-6' + 30)], where the chirality 
and vorticity are parallel. Here, A±{r) are real, A_|_(r oo) = Abcs(2^)) and 
A_(r ^ oo) = 0. 

We focus on SQ{R — 0), where both the defect and the vortex center are 
situated just at the origin r = 0. The quasiparticles inside the vortex core, 
subject to the Andreev reflection, run along straight lines, namely the quasi- 
classical trajectories [11,12]. We consider the quasiclassical trajectories which 
go through the origin r = 0. On those trajectories with zero impact param- 
eter, the position vector is parallel to the direction of the trajectory (i.e., 
r II k), and therefore (j) = 6, ^ -|- tt. In this situation the order parameter 
is A-^- {r, (p; 9) = ±(A+(r) + A_{r)) in the case of (i), and A^+(r,0;6') = 
±(A_|_(r) -|- A_(r)) exp(2i6') in the case of (ii). Of importance is the difference 
in the phase factor of these order parameters of (i) and (ii). 

On the basis of an analysis of the so-called zero-core vortex model in Ref. [3], 
the matrix elements of ^imt at the vortex center are approximately obtained 

as Qimt = \JujI + |A|Va;„, /i^t = -A/un, and f^^^ = A*/uJn, where A = 
Aj|"^(r — >■ oo,0 = 9;9). Because of the difference in the phase factor of the 
order parameters mentioned above, (/imt)6> — /imt and {f}^i)e = /iLt the 
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case of (i), while (/imt)6» = and (/imt)e = in the case of (ii). Therefore, in 
the case of (i), (^imt)6» = ^imtj and then [t,^imt] — from Eq. (5). In the case 
of (ii), (^imt)^ 5imt, and then [i,gi^t] 0. 

When [t 

, 9imt] — 0, the Eilenberger equation (4) in the presence of the defect is 
identical to Eq. (2) (the equation in the absence of the defect), i.e., the defect 
has no influence on the Green function and the free energy. From it and the 
above analysis, 5^2(0) = in the case of (i) when the chirality is antiparallel to 
the vorticity, and 6il{0) 7^ in the case of (ii) when the sense of the chirality 
is the same as that of the vorticity. It means that the vortex pinning depends 
on the chirality in the chiral p-wave superconductor. This analytical result 
is certainly conflrmed by a numerical calculation [13] with self-consistently 
obtained order parameters A^^(r, 0;6') [11]. 



4 Summciry 

We investigated the elementary vortex pinning potential 6Q{R) on the basis 
of the quasiclassical theory of superconductivity. The numerical results were 
presented in Fig. 1 for the d-wave and s-wave pairings. The vortex pinning 
energy in the d-wave pairing case is about 10 times larger than those in the 
s-wave pairing cases, at a high temperature, because of the breakdown of the 
Anderson's theorem in unconventional superconductors. We also investigated 
analytically the quasiparticles inside the vortex core in the chiral p-wave su- 
perconductor, and found the chirality-dependent vortex pinning. We expect 
that the chirality-dependent vortex pinning revealed in this paper shares the 
same physics with impurity effects on the quasiparticle relaxation time [12] 
and on the energy spectrum [14,15] inside the vortex core in chiral p-wave 
superconductors. 
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